Coherent transport in a two-electron quantum dot molecule 
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We investigate the dynamics of two interacting electrons confined to a pair of coupled quantum dots 
driven by an external AC field. By numerically integrating the two-electron Schrodinger equation in 
time, we find that for certain values of the strength and frequency of the AC field the electrons can 
become localised within the same dot, in spite of the Coulomb repulsion between them. Reducing 
the system to an effective two-site model of Hubbard type and applying Floquet theory leads to 
a detailed understanding of this effect. This demonstrates the possibility of using appropriate AC 
fields to manipulate entangled states in mesoscopic devices on extremely short timescales, which is 
an essential component of practical schemes for quantum information processing. 
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A set of electrons held in a semiconductor quantum 
dot is conceptually similar to a set of atomic electrons 
bound to a nucleus, and for this reason these structures 
are sometimes termed "artificial atoms" [[[]. To extend 
the atomic analogy further, we can consider joining quan- 
tum dots together to form "artificial molecules". The 
simplest example consists of two coupled dots containing 
a single electron. The dynamical properties of this sys- 
tem when driven by an AC field have been extensively 
studied (|] , and the technique of Floquet analysis || nas 
proved to be particularly effective. The pioneering work 
in Ref. Q first noted the dramatic effects on the tunnel- 
ing rate arising from crossings and anti-crossings of the 
Floquet quasi-energies, and the term dynamical localisa- 
tion was coined to describe the phenomenon in which the 
AC field appears to trap the electron within just one dot 
over long time scales. Subsequent analytic and numerical 
work ||, Q showed that for high frequencies this locali- 
sation occurs when the ratio of the field strength to the 
frequency is a root of the Bessel function Jo. 

Adding a second electron to the coupled dot system 
introduces considerable complications, as for realistic de- 
vices the Coulomb interaction between the electrons can- 
not be neglected, and a single-particle analysis is not ap- 
plicable. Achieving an understanding of the dynamics of 
this system is, however, extremely desirable, as the abil- 
ity to rapidly control the localisation of electrons using 
AC fields immediately suggests possible applications to 
quantum metrology and quantum information process- 
ing. In particular, manipulating a pair of entangled elec- 
trons on short timescales is of great importance in the 
rapidly developing field of quantum computation. Al- 
though numerical investigations have been made of the 
rich behaviour displayed by interacting electron systems 
driven by an AC field 0, there is, at present, little under- 
standing of the observed effects beyond phenomenology 
and numerical experimentation. 

We address this problem here by applying the Floquet 
formalism to a system of interacting particles. We ini- 
tially consider a realistic model of two interacting elec- 
trons confined to a pair of coherently coupled quantum 



dots, and study its response to AC fields by a numeri- 
cal method in which the Coulomb interaction is treated 
exactly. We find the surprising result that even in the 
presence of strong inter-electron repulsion, suitable AC 
fields can localise both electrons within the same dot. 
We further show that this effect is also exhibited by a 
simple two-site model, obtained as a reduced form of the 
quantum dot system. As in the case of the single-particle 
model, we find Floquet theory to be an extremely power- 
ful tool, and we succeed in generalising the single-particle 
solution ||, to include interactions, from which we can 
easily and accurately identify the parameters of the AC 
field for which localisation occurs. 

We study a quasi-lD structure, in which the linear 
dimensions of the quantum dots in the (x, y) plane are 
much smaller than their length. The energies associated 
with transverse excitations are thus much higher than for 
longitudinal motion, and can be neglected. The dots are 
separated by a potential barrier which regulates the de- 
gree of tunneling between the dots. Strongly coupling the 
dots allows the electrons to form extended states ( "molec- 
ular orbitals") over the entire double-dot system, which 
can be explicitly seen in transport experiments ||, ||. 
Within the effective mass approximation, the Hamilto- 
nian for the system is given by: 
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where z\, z 2 are the spatial coordinates of the two elec- 
trons, m* is the effective mass and e is the effective 
permitivity. GaAs material parameters are used, with 
to* = 0.067to c and e = 10.9e - E(t) is the external AC 
field, E(t) — Ecos(ujt). To simplify the numerical treat- 
ment, the system was placed within an infinite square 
well of length L to prevent the electronic wave-function 
from leaking out of the dot structure. So that these in- 
finite barriers did not unduly influence the properties of 
the quantum dots, thick buffer zones were included be- 
tween these barriers and the walls of the quantum dots 
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so that the electronic wavefunction effectively decays to 
zero before reaching the barriers. The confining poten- 
tial, V(z), is plotted in Fig.|l|. The Coulomb potential 
used, Vc{r), had the form: 

where A is a measure of the transverse width of the struc- 
ture (taken to be 1 nm in this investigation). This form 
for Vc reproduces the normal 1/r fall-off of the Coulomb 
potential, but is non-singular in the limit r — ► 0. 

To study the time evolution of the system we used the 
eigenstates of the static Hamiltonian, H , as a basis, since 
expanding in these states is a well-controlled procedure, 
and also ensures that correlation effects arising from the 
Coulomb interaction are automatically encoded within 
each basis function. We note that as Eq.[l] contains no 
spin- flip terms, there is no mixing between the singlet and 
triplet sub-spaces. In particular, if the initial state has a 
definite parity this symmetry is retained throughout its 
time evolution, and consequently only basis functions of 
the same symmetry need to be included in the expansion. 
To obtain the eigensystem of either the singlet or triplet 
sub-spaces, we employed the Lanczos technique described 
in Ref. [(To), which is particularly efficient for systems in 
which the interaction is diagonal in real-space. 

The initial state used was the ground state of H (a 
singlet). In the absence of the external electric field this 
would have a trivial time evolution, simply acquiring 
a phase. Applying the field, however, causes the ini- 
tial state to evolve into a superposition of eigenstates: 
IV'M) = J2 c n(t) \E n ). Substituting this expansion into 
the time-dependent Schrodinger equation yields a first- 
order differential equation for the expansion coefficients 

Cn ■ 

dc n M 

fa— = c n {t)E n - eE(t) 22 F nm c m (t), (3) 

m— 1 

where E n is the n-th eigenvalue of H , and F mn are 
the overlap integrals of the dipole operator, F mn = 
(E m \(zi + z 2 )\E n ) . A fourth order Runge-Kutta method 
was used to time-evolve the expansion coefficients using 
Eq|| The number of basis states used in the expansion, 
M, required for convergence depended on the strength of 
the electric field, and values of M = 50 were required for 
the strongest field considered. 

A similar model was studied recently in Ref. for a 
larger system size of L = 100 nm. An important ad- 
vantage of considering a smaller system, however, is that 
finer structure can be resolved in the results due to the 
larger spacing between the energy levels. Following Ref. 

we define a conditional probability function, which 
gives an extremely useful description of the state of the 
system. The probability that one particle is in the right 
dot while the other is in the left is given by P RL (t): 



P RL (t) = 2 f dz x [ dz 2 \^(z u z 2 )\ 2 , (4) 
Jr Jl 

where the notation R / L signifies that the integration 
is taken over the right/left quantum dot. P RL takes a 
value of one for maximally delocalised states (when one 
electron is in the right dot and the other is in the left), 
and is zero for localised states when both electrons oc- 
cupy the same dot. The initial state plotted in Fig.|l| 
is highly delocalised, and has the value P RL = 0.849. 
Our investigation consists of evolving this state through 
a time period of 18 ps while measuring P RL (t). We term 
the minimum value of P RL attained during this period 
-Pmin ■ and use this to quantify the degree to which the 
AC field brings about localisation. 

In Fig.^ja we present a contour plot of P m in as a func- 
tion of the parameters of the AC field. Dark areas cor- 
respond to low values of P m i n , indicating a high prob- 
ability that both electrons are occupying the same dot. 
Surprisingly this can occur even at weak field strengths, 
despite the presence of the Coulomb interaction. We 
note that the dark areas form horizontal bands, indicat- 
ing that for various "resonant" values of u>, localisation 
can be produced over a wide range of fields. The spacing 
of the bands decreases with u>, and for values of hcu < 2.8 
meV the structure is too fine to be resolved. Between 
these bands strong localisation is not produced. These 
results are qualitatively similar to those of Ref. R] but 
show finer detail, allowing us to additionally observe that 
these bands are punctuated by narrow zones in which the 
field does not create localisation. Their form can be seen 
more clearly in the cross-section of P m i n given in FigJ^a, 
which reveals them to be narrow peaks. These peaks are 
approximately equally spaced along each resonance, the 
spacing increasing with oj. 

We emphasise that these results are radically different 
to those obtained for non-interacting particles. In this 
case an analogous plot of delocalisation shows a fan-like 
structure Q , jllj , in which localisation occurs along lines 
given by to = E/xj , where Xj is the j-th root of the Bessel 
function Jq(x). As a test of our method, we repeated 
our investigation for a quantum dot system in which the 
inter-electron Coulomb repulsion was set to zero, and 
found that the fan structure was indeed reproduced. 

To account for these results it is thus necessary to go 
beyond the non-interacting case. We consider a highly 
simplified model in which each quantum dot is replaced 
by a single site. Electrons can tunnel between the sites, 
and importantly, we include interactions by means of a 
Hubbard ?7-term: 
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H = -tJ2 ( c Lc 2CT + h.c.) + {Un^nn + E^tjm) . 
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(5) 

Here t is the hopping parameter. In this analysis we set 
% = 1 and measure all energies in units of t. Ei(t) is the 
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external electric potential applied to site i. Clearly only 
the potential difference, E\ — E2, is of importance, so we 
may choose the convenient parameterisation: 



Exit) = — COSLUt, 



E 2 (t) = --cosujt. (6) 



A numerical investigation of such a model was made re- 
cently in Ref. |l2j for the case of very weak fields. The 
Hilbert space of Hamiltonian ([5]) is six-dimensional, com- 
prising three singlet states and a triplet. As with the 
double-dot system, the singlet and triplet sub-spaces are 
completely decoupled, allowing us to study just the sin- 
glet space. In the absence of the AC field the eigenvalues 
of the singlet Hamiltonian can be found analytically, and 
for large U they consist of two almost degenerate excited 
states, separated from the ground state by the Hubbard 
gap U. This mimics the eigenvalue structure of the low- 
est multiplet of states of the full double-dot system. We 
time-evolve the system following the same procedure as 
before, using the ground state as the initial state. As 
the singlet Hamiltonian is only three-dimensional, how- 
ever, the computations can be done correspondingly more 
rapidly. In Fig.|]b we show the results for P m i n obtained 
by setting U — 8, which clearly shows that this simpli- 
fied model strikingly reproduces the behaviour of the full 
system of Eq.[l]. 

As the driving field (|^) is a periodic function of time, 
we can make use of Floquet analysis to describe the time 
evolution of the system in terms of its Floquet states and 
quasi-energies ||. The Hamiltonian (||) is invariant un- 
der the combined parity operation x — ► — X] t — ► t + T/2, 
and so the Floquet states can also be classified into these 
parity classes. Close approaches of the Floquet quasi- 
energies as the system parameters are varied produce 
large modifications of the tunneling rate, and hence of 
the system's dynamics. Quasi-energies of different parity 
classes may cross, but if they are of the same class they 
form an anti-crossing. Numerically the quasi-energies 
can be conveniently obtained by diagonalising the time 
evolution operator for one period of the driving field, 
U(t + T, t). This is particularly suited to the numeri- 
cal approach we have used, as U(t + T,t) is simply the 
operator obtained by time-evolving the identity matrix 
J3 over one period of the driving field. 

We present in Fig.||a the Floquet quasi-energies as a 
function of the field strength for oj = 2, one of the reso- 
nant frequencies visible in Figj^b. Below the spectrum we 
also plot the behaviour of P m i n . We see that the system 
possesses two distinct regimes of behaviour. For weak 
fields, as studied in Ref. |L2]], the Floquet spectrum con- 
sists of one isolated state (which evolves from the ground 
state) and two states which make a set of exact cross- 
ings. In this regime -Pmin decays slowly to zero, showing 
little structure. As the field strength exceeds U, however, 
this abruptly changes to a novel, previously unseen, be- 
haviour in which P m i n remains close to zero except at 



a series of narrow peaks, corresponding to the close ap- 
proaches of two of the quasi-energies. A detailed exam- 
ination of these approaches (see Fig.||b) reveals them to 
be anti- crossings between the Floquet states which evolve 
from the ground state and the higher excited state, and 
have the same parity. The remaining state, of opposite 
parity, makes small oscillations around zero, but its exact 
crossings with the other two states do not correlate with 
any structure in P m i n . 

To interpret this behaviour we seek analytic expres- 
sions for the quasi-energies. We choose to use a perturba- 
tional approach S, starting from the Floquet equation: 



H 



d 



(7) 



Here H is the full Hamiltonian (||), ej are the Floquet 
quasi-energies, and \4>j(t)) are the Floquet states. Our 
procedure is to first find the eigenstates of the operator 
(Hi — where Hi is the second term in EqJ|| con- 

taining all the interaction terms, and then treat the tun- 
neling component H t as a perturbation. An important 
advantage of this approach is that the Floquet states are 
stationary states of (0). Consequently, by working in an 
extended Hilbert space of T-periodic functions |l3| , the 
corrections can be evaluated easily by standard Raylcigh- 
Schrodinger perturbation theory, without requiring more 
complicated time-dependent methods. 

In a real-space representation the interaction terms are 
diagonal, and so it can be readily shown that an orthonor- 
mal set of eigenvectors of (Hj — ij^) is given by: 



|e (i)) = (exp[ieot], 0, 0) 



\e+{t)) = (0, exp 



E 



-i(U — e+)f + i — sinujt 



0) 



|e_(i)) = (0, 0, exp 



E 

-i(U — e-)t — i — smut 
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Imposing T-periodic boundary conditions reveals the cor- 
responding eigenvalues (modulo lu) to be e = and 
e± = U. These eigenvalues represent the zeroth-order 
approximation to the Floquet quasi-energies, and for fre- 
quencies such that U — n lo all three eigenvalues are de- 
generate. This degeneracy is lifted by the perturbation 
H t , and to first-order the quasi-energies are obtained by 
diagonalising the perturbing operator Py = ((ei\H t \e.j)) , 
where ((. . .)) denotes the inner product in the extended 
Hilbert space || . By using the identity: 

00 

exp [— i/3 sin ujt] = J m {P) exp [— imujt] , (9) 

m— — 00 

to rewrite the form of |e±(t)), the matrix elements of 
P can be obtained straightforwardly, and its eigenvalues 
subsequently found to be eo = and e± = ±2J„(E/uj). 
This solution clearly reduces to the well-known solution 
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for non- interacting particles when {7 = 0. Fig.|3|a demon- 
strates the excellent agreement between this result (with 
n = 4) and the exact quasi-energies for strong and mod- 
erate fields, which allows the position of the peaks in 
P m i n to be found by locating the roots of J n . Similar 
excellent agreement occurs at the other resonances. For 
weak fields, however, the interaction terms do not dom- 
inate the tunneling terms and the perturbation theory 
breaks down, corresponding to the weak-field regime in 
which P m i n decays smoothly to zero. Away from the res- 
onances the first-order correction identically vanishes, re- 
sulting in the lack of structure seen between the resonant 
bands in Fig.|]b, and it is necessary to go to higher or- 
ders in perturbation theory to obtain the quasi-energy 
behaviour. 

In summary, we have investigated the dynamics of an 
interacting two-electron system driven by an AC field. 
We find that despite the Coulomb interaction a suitable 
AC field can nonetheless produce localised states. This 
localisation occurs over a range of field strengths at fre- 
quencies for which an integer number of quanta, n, is 
equal to the interaction energy. At high fields we find 
a novel regime of behaviour in which this localisation 
vanishes at the roots of J n (E/w), and explain this using 
perturbation theory. These results are of general appli- 
cability to AC-driven systems of interacting electrons, 
and hold out the exciting prospect of controlling and 
manipulating correlated quantum states on picosecond 
timescales by means of applying AC fields. 
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FIG. 1. Geometry of the double-dot system. The dotted 
line plots the charge density of the ground state wavefunction. 
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FIG. 2. -P m ; n as a function of the strength E and energy Tioj 
of the AC field: (a) for the interacting quantum dot system 
(Tkj in units of meV) (b) for the Hubbard model with (7 = 8 
(both axes in units of t). 
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FIG. 3. (a) Quasi-energy spectrum for the two-site model 
for (7 = 8 and uj — 2, circles = exact results, lines = pertur- 
bation theory, (b) magnified view of exact results for a single 
anti-crossing. Beneath are the corresponding plots of P m ; n . 
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